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A new representation for electrons is introduced, in which the electron operators are written in 
terms of a spinless fermion and the Pauli operators. This representation is canonical, invertible 
and constraint-free. Importantly, it simplifies the Hubbard interaction. On a bipartite lattice, 
the Hubbard model is reduced to a form in which the exchange interaction emerges simply by 
decoupling the Pauli subsystem from the spinless fermion bath. This exchange correctly reproduces 
the large U superexchange. Also derived, for U = ±oo, is the Hamiltonian to study Nagaoka 
ferromagnetism. In this representation, the infinite-?/ Hubbard problem becomes elegant and easier 
to handle. Interestingly, the ferromagnetism in Hubbard model is found to be related to the gauge 
invariance of the spinless fermions. Generalization of this representation for the multicomponent 
fermions, a new representation for bosons, the notion of a 'soft-core' fermion, and some interesting 
unitary transformations are introduced and discussed in the appendices. 



PACS numbers: 71.10.Fd,75.10.Jm,05.30.Fk 



I. INTRODUCTION 

Strongly correlated electrons represent a large class 
of physical systems exhibiting a variety of novel phe- 
nomena 1 ' 2 . The transition metal oxides, such as the 
higli-Tc superconducting cuprates 3,4,5,6 , or the colos- 
sal magneto-resistive manganites 7 , provide some of the 
most famous and actively investigated materials of this 
kind. The common feature in all such materials is the 
dominant presence of local repulsion between electrons, 
with manifestly non-trivial consequences, such as the 
antifcrromagnctic-insulating behavior at half-filling (the 
Mott insulators), which can not be understood from the 
band-theory of solids 8,9 . Physics of these systems is in- 
trinsically non-perturbative, and hence difficult. The 
Hubbard model 1,10,11 ' 12,13 , which describes a system of 
tight-binding electrons with on-site repulsion, allows the 
simplest meaningful description of the strongly corre- 
lated electrons. Since the local repulsion U (also called 
the Hubbard interaction) and the hopping t (assuming 
nearest- neighbor) have opposing effects on the electrons, 
the ground state of the Hubbard model can be metal- 
lic/insulating, magnetic/non-magnetic or something else 
(yet unexplored), depending upon the ratio t/U, the elec- 
tron density n e , and the geometry of the lattice. 

For sufficiently large U at half-filling (n e = I), the 
Hubbard model is known to systematically produce the 
antifcrromagnetic exchange interaction and the insulat- 
ing ground state 8 . As t/U is increased, the insulat- 
ing phase undergoes a quantum transition to a metallic 
phase. Though there have been many elaborate studies 
on Hubbard model 9,14 , there is no final verdict as yet on 
the nature of this metal-insulator transition (except in 
one dimension where the exact ground state is insulat- 
ing for any non-zero U 15 ). Also, the search for a super- 
conducting ground state in purely repulsive lattice mod- 
els 16 ' 17,18 has become very important since the discovery 
of superconductivity in cuprates and other strongly cor- 
related materials, but with no definite answer. Present 



understanding of the Hubbard model is far from being 
complete, and has been a subject of great research inter- 
est 13 . The main difficulty lies in the handling of the Hub- 
bard interaction. This has led to the development of nu- 
merous approximate methods. The Gutzwillcr projected 
variational wavefunctions 19,20 , the dynamical mean-field 
theory 14 , the slave-boson (fermion) approach 21,22 or the 
cluster approximations 23 are a few important ones to 
mention. Exact numerical diagonalization has also been 
used to study Hubbard model on small lattices 24 . The 
motivation to search for newer and more effective schemes 
of investigation, for a better understanding of the Hub- 
bard model (and the strongly correlated systems in gen- 
eral), is therefore very strong. 

In this paper, we present a new canonical represen- 
tation for electrons, and demonstrate its usefulness for 
the Hubbard model. The basic idea is to identify the 
independent physical attributes of an electron, in terms 
of which the electron can be described as a composite. 
The two attributes of an electron are its spin and the 
Fermi-ness. The spin defines the two-level structure of 
the electron which can be described using Pauli oper- 
ators. The Fermi attribute, which assigns the correct 
quantum statistics (anticommutation) to the electron, is 
described in terms of a spinless fermion. In our descrip- 
tion, these two canonical degrees of freedom are com- 
bined to make the electron. An important feature of this 
representation is that it transforms the Hubbard inter- 
action (the particle-hole symmetric form) into the chem- 
ical potential for the spinless fermions. The simplifica- 
tion thus achieved further encourages us to investigate 
the Hubbard model in this representation. We discuss 
two very important cases, that is the superexchange for 
finite U and the Nagaoka ferromagnetism 25 for infinite 
U, to highlight the physical usefulness of our represen- 
tation in describing strongly correlated electrons. The 
paper is organized as follows. In the next section, we in- 
troduce the representation, discuss its canonical nature, 
and derive the inverse representation. In section III, we 
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TABLE I: Local Hilbert spaces and the operators of the elec- 
tron, spinless fermion and a Pauli (two-level) system. Here, 
|0) and |0) denote the vacua of the electron and the spinless 
fermion, respectively. Operators f a and & satisfy the canoni- 
cal Fermi algebra, while cr z and <x are the Pauli operators. 
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discuss the Hubbard model. On a bipartite lattice, we 
reduce the model to a nice form in which a part of the 
hopping explicitly generates the exchange interaction in a 
simple decoupling scheme. For large U at half-filling, this 
exchange interaction correctly reproduces Anderson's su- 
perexchange 8 . In section IV, we present the Hamiltonian 
formulation of the Nagaoka ferromagnctism. In the ap- 
pendices, we discuss a few related but distinct topics. 
First, we generalize this representation for A^-component 
fermions in appendix A. Next, we write down a new rep- 
resentation for the bosons. In appendix C, we introduce 
a unitary transformation with Majorana gauge, and de- 
rive the Jordan- Wigncr transformation 26 ' 27 . Finally, we 
introduce the notion of a 'soft-core fermion', and discuss 
its basic properties. 



II. NEW REPRESENTATION FOR 
ELECTRONS 

Consider the Hilbert space of electron on a single site. 
The four basis states which span this local electronic 
Hilbert space are tabulated in Table I, where |0) denotes 
the empty site, and /j and /j are the creation operators 
of T and i electrons, respectively. Also consider a spin- 
less fermion and a two-level (Pauli) system, each hav- 
ing a two-dimensional Hilbert space (see Table I). Let a 
and be the annihilation and creation operators of the 
spinless fermion, and a z and er^ be the usual Pauli opera- 
tors describing the two-level system. Since the composite 



TABLE II: A mapping between the states of the electron and 
a composite system of a spinless fermion and a Pauli system. 
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Hilbert space of a spinless fermion and a Pauli system is 
four-dimensional, we can define a one-to-one mapping be- 
tween the basis states of this composite Hilbert space and 
that of the electron 28,29 . In Table II, we define one such 
mapping, in which the one-electron states correspond to 
having one spinless fermion, and the empty or double- 
occupancy electronic states correspond to the vacuum 
of the spinless fermion. For the one-electron states, |+) 
corresponds to f spin, and |— ) to J, spin of the electron. 
Otherwise, they correspond respectively to the double- 
occupancy or empty state. This mapping, which is mo- 
tivated by our idea of the electron as a composite of the 
Fermi and Pauli degrees of freedom, results in the follow- 
ing new representation for the electron operators. 

n = 4><? + (i) 



Here, 4> = + a, and i'tp — a* —a. The operators (f> and 
(called the Majorana or real fermions) arc Hcrmitian, and 
satisfy the relations: <f> 2 = if 2 = 1 and {<f>, ip} = 0, where 
{, } denotes the anti-commutator. As usual, the spinless 
fermion operators commute with the Pauli operators. We 
have derived Eqs. (1) and (2) using the definition of fl 

and f} in terms of the basis states (as given in Table I). 

In this representation, the electron number operators 
are given as: n-j- = (1 + c z )/2 and nj, = | + — n) a z , 
where h = ala. The physical spin, S, of the electron, 
which is defined as S z = (n-f — fi\)/2 and S + = ftfi, 
takes the form: S = na/2. Similarly, the local pairing 
(the pseudo-spin operators), P can be written as: P = 
(n - l)oy2, where P+ = /|/| and P z = (n T + n l - l)/2. 
Also, S + P = a/2 and SiPj = 0, where i, j = x,y,z. 
Here, a x = a + + a~ and a v = —i{a + — <j~). Physi- 
cally speaking, a represents the spin of electron in the 
presence of the spinless fermion. Otherwise, it represents 
local pairing. These forms (especially that of S) are con- 
sistent with our intuitive picture in which an electron is 
envisaged as a "composite" object made up of the Fermi 
and the spin attributes. Furthermore, it is important 
to note that (n^ — |) — 5) = 5 (3 — n), that is the 
Hubbard interaction between electrons transforms into 
the number operator of the spinless fermion. 

The new representation of the electron is canonical, 
constraint-free and invcrtible. It is canonical because the 
electron operators satisfy the usual anti-commutation re- 
lations. It is obviously constraint-free because there are 
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no extra conditions to be satisfied by a or a operators. 
This is a consequence of the one-to-one mapping between 
the two Hilbert spaces. Since there are no excess or lost 
states, there is no constraint. Also, the spinless fermion 
and the Pauli operators can be represented uniquely in 
terms of the electron operators. The inverse of the new 
representation is given as: 

« f = (1 (3) 
g z = 2n T - 1 (4) 

* + =f](f! + h) (5) 

The inverse representation is also consistent, because the 
fermion operator, a, satisfies the anti-commutation rela- 
tion and the Pauli operators satisfy the angular momen- 
tum algebra. The Pauli operators also satisfy the anti- 
commutation relation: {<7 4 ,(7 J } = 25y for i,j = x,y,z. 
And of course, a and a' commute with a. 

The inverse representation gives us more insight into 
the nature of the new objects. Equations (4) and (5) 
are same as the Majorana fermion representation of the 
spin- 1/2 operators 30 . The spinless fermion [Eq. (3)] is ef- 
fectively a Gutzwillcr-projccted electron whose creation 
is equivalent to the creation or annihilation of j elec- 
tron, subject to the absence or presence of f electron 
respectively. It is also like the Bogoliubov quasiparticle 
in the BCS theory as it mixes a particle with a hole. 
Surprisingly, the inverse representation has more famil- 
iar appearance, and hence, it is easier to appreciate. We 
can generate many equivalent forms of our representa- 
tion by making suitable unitary rotations. This gives us 
the freedom of choosing a form that may serve us better 
for a given problem. For example, the unitary operator, 
U = 1 — 2h'[hi, transforms Eqs. (3-5) as follows. 

U&W = St = (l-n T )/t + n T / i (6) 
Ua+U^ = r+ = -/f(/[ _ /J ( 7 ) 

Or conversely, fj — iij>bT + , and /j = (<fo — i^if, t z )/2, 

where ^ — w + b and if/ty = w — b. We will use this and 
the original form of the representation to conveniently 
describe the Hubbard model on a bipartite lattice. 

III. HUBBARD MODEL AT HALF-FILLING 

Consider the Hubbard model with nearest-neighbor 
hopping on a lattice. Its Hamiltonian in terms of the 
electron operators can be written as: 

H = ~ 2 (/r,s/r+£,s + /r+^s/r.s) 

r,S s=T,l 




where r is summed over all the sites of a lattice, and 
5 is summed over the nearest-neighbors of a site. The 



electron hopping is denoted by t, and U is the Hubbard 
interaction. In Eq. (8), the factor of 1/2 in t is there 
to undo the double-counting of the hopping due to the 
explicit use of its Hermitian conjugate. Therefore, the t 
we use is same as the one commonly reported. The U- 
term in Eq. (8) is written in the particle-hole symmetric 
form 31 . The Hubbard model in the new representation 
[using Eqs. (1) and (2)] takes the following form. 

H = ~ J X] ®r &r ~ \ a r+S - h - C ) 4>r ^+8 

~\ Yl + CT* + 5 k+5_ Vv) + jL (9) 

r,S 

Here, L is the total number of lattice sites. In this form, 
the Hubbard interaction has taken the role of chemical 
potential for the spinless fcrmions, and the hopping has 
transformed into a complicated spin-fcrmion interaction. 
The Hubbard model has thus turned into a problem of 
spin-1/2 objects immersed in a grand-canonical bath of 
spinless fermions, which are constantly being created and 
annihilated around an average fermion density decided by 
U. The total a z is however conserved, which is same the 
conservation of the total number of "["-electrons. 

Since the local gauge of the electron governs the flow 
of charge-current in a system, it is important to consider 
how the new objects respond to such a local gauge trans- 
formation. This will give us an idea of the way charge 
is carried by the new objects. Let /J s — > e 1 r /J a be 
the gauge transformation on electrons, under which the 
spinless fermions and the Pauli operators transform as: 

at - (l-n rT )/ r V^-n rT / ri e^ (10) 

-r + - /rV^ + A/rl (H) 

This implies that a spinless fermion carries one unit of 
electron or hole charge (of J, electron, to be more specific) , 
and <7+ carries a charge of two electrons or no charge at 
all. It can be understood by explicitly considering the ac- 
tion of and a + on the electronic states. For complete- 
ness, we also write the current density operator for the 
Hubbard model in terms of the spinless fermions and the 
Pauli operators, which is J T = (iet/2) J2s §.(fr,sfr+5 s ~ 

fr+S,Jr,s) = (i^/ i )T,l^.bl ) Ti>T+S + 4>r<t>r+s(Vr ' &r+s)]- 

Interestingly, the current operator has nicer form as com- 
pared to hopping in the new representation, and it also 
reveals the exchange interaction. At this point, it may 
be mentioned that Eq. (6) is same as Eq. (10) (with some 
fine tuning) for 9 r = ir/2. In fact, the U in Eq. (6) was 
suggested by this observation. 

On a bipartite lattice, Eq. (9) can be transformed to a 
more elegant form by making a suitable unitary rotation 
on the operators on one sublattice. Let A and B denote 
the two sublattices. We write the electron operators on 
A-sublattice, in terms of a, and <?, as /| = a c + 
and /| = (iipa — 4>aC z )/2. The electrons on B-sublattice 
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are described in terms of b, and f such that ft = 

itpb T+ an d f \ = {4>b ~ i^& rZ )/2- The representation on 
A-sublattice is in the original form [as in Eqs. (1) and (2)], 
while on B-sublattice, it is a unitary-transformed version 
of the same [see Eqs. (6) and (7)]. The Hubbard model 
on a bipartite lattice can therefore be written as: 



H = 



—L — — 

4 2 

it 



Y a R a R + Y ^r 6 r 



ReB 



(12) 



^£ r 

2 Y Y [^a^b.R+S + ^b,R+s4>a,R (°R ' ^R+i) 
ReA 6 

This form of H is neat and suggestive 32 . It presents 
the dual role of hopping very clearly. While a part of 
the hopping contributes purely to the dynamics of the 
spinlcss fcrmion bath, the remaining part generates an 
SU(2) invariant exchange interaction 33 , whose strength 
is determined by the fluctuations in the bath. 

To gain further insight, we discuss the Hubbard model 
in a decoupled form: H — H Fe 

U 



H PauW whCTC 



Y "r" r + T Y b ^R 



(13) 



ReA 



ReB 



~ f YY ( a R ^R+5 + ^R+S «R 
R£A 8 

u ■ 



~ tz Y (jk "ik^-k + 7k ^-kflk 



(14) 
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FIG. 1: Antiferromagnetic exchange J as a function of U. The 
superexchange interaction, 4t 2 /U (t = 1 in our calculation), 
is also plotted for comparison. For large U, all J — » 4t 2 /U. 



z S<5 



Jk-S 



7k | e k , and z is the nearest-neighbor co- 



ordination. Since 7^ = 7-k, therefore, |7_k| = l7k| and 
A_k = — ^k- Here, the wavevector k lies in the sublat- 
tice Brillouin zone, and the Fourier transformation of the 



spinlcss fcrmions is defined as: 



E'zk-R ' 
k e °k, 



l R — Y L ^k 

and similarly for 6k- The Hpermi is diagonalized by ap 
plying two unitary rotations in following way. 



Fermi — 



- 7/ 7/ 



Ug U x H f 
y^^k (1 - a k «k - b"\J> k 



(16) 



and 



H 



Pauli 



K 

2 



YY ( ff B. ■ TR+8 + 1 



ReA s 



(15) 



This decoupled from is derived in the following way. 
First, we rewrite er R • f R+5 in Eq. (12) as —1 + (1 + cr R • 

r R+(5 ), and include —ip b R+( 50 a R resulting from this into 
the purely spinlcss fcrmion part of H , which is now called 
H Fermi- Then, replace ^bR+s&aR [which is in multi- 
plication with (1 + cj r • Tr, 5 )] by its expectation value 
in the ground state of H Fermi- The resultant involving 
only Pauli operators is called Hp a uU- Thus, the bath of 
spinlcss fcrmions in the decoupled picture is described by 
H Fermi, and H Pauli is the exchange Hamiltonian, where 
(<?r • TR+s + l)/2 is the Dirac-Hcisenberg exchange oper- 
ator on a nearest-neighbor bond. Note that the constant 
term UL/A has been absorbed in rewriting 6 R S R as 6 R 6 R . 

The validity of this decoupling can be tested against 
the well-known superexchange interaction, At 2 /U, in the 
limit of large U. The strength of the exchange interac- 
tion, J, is given as: J = 2t(, where ( = -i{i>b,R+54>a,R) 
(uniform average on each nearest-neighbor bond) is de- 
termined purely by Hpermi- Equation (14) is the mo- 
mentum space representation of Hpermi, where 7k = 



The unitary operators in Eq. (16) are defined as: U\ 



H k e ix ^ bk , and Ug = J~I k e tiki ^ b '-^- b - khk \ where U x 
absorbs the phase Ak into 6k , and Ug diagonalizes the 
dk, &L k mixing. The quasi-particle dispersion, £ k = 

1 + (z£|7k|) 2 , is positive, and tan 20k = zt|7k|/(y)- 
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The ground state of Hpermi is given by d k d 
for all k. We calculate ( in the ground state of Hp 



and find that J 



~£~ Sk ' which is certainly pos- 
itive. Hence, the exchange interaction in Hp au n is anti- 
ferromagnetic. Explicitly on a £)-dimensional hypercubic 
lattice, where z = 2D and 7k = 52i=i cosfc^, we cor- 
rectly retrieve Anderson's superexchange, J — > 4i 2 /|[/|, 
for U 2zt. In Fig. 1, the calculated values of J are 
plotted as a function of U/t, for a few different bipartite 
lattices. The densities of a and b type spinless fcrmions 
are calculated to be: n a = Tib = n = ^ + For 
U > 2zt, n — > 1, and for U — > 0, n — ► | (see Fig. 2). 
Both of these limiting cases give the expected values of 
the spinless fermion density. Since the local densities 
n a ,R and nb iR also take the uniform value n, the electron 
density can be calculated as: n e = 1 + ^ r cr^, which 



is equal to 1 when total a z = (that is, J2a crZ +Ss 7-2 = 
0). Therefore, in the decoupled problem, the half-filling 
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FIG. 2: Density of the spinless fermions, n a — rib — n vs - U. 
For U 3> 2zt, n ~ 1, and it is 0.5 for U = (which is the 
correct value for the non- interacting case). 



for electrons is achieved when the number of f electrons 
is exactly L/2. Since the ground state of Hp a uii cor- 
responds to total a z = 0, the half-filled condition for 
electrons is always satisfied in the ground state of the 
decoupled problem. Hence, the decoupling scheme dis- 
cussed here is a consistent approximation, as it correctly 
yields the supcrcxchangc model of electron spins at half- 
filling for large U (S = ha/2 ~ a/2 for (n) ~ 1). 

We close this section with the following observations. 
It is perhaps the first time that the supcrexchange in- 
teraction (which is truly a strong correlation effect) has 
been derived in a mean-field like approximation. This 
is to be contrasted with the usually known approaches 
which give rise to the supcrexchange only in the second 
order perturbation theory 8,22 . Another point to be noted 
is that J is an even function of U. Therefore, the electron 
pseudo-spins [P = (1 — h)a/2 ~ a/2 for (h) — > 0] also 
interact via antiferromagnetic 'superexchange' for large 
negative U. This can be deduced, without explicit cal- 
culation, from the particle-hole (ph) transformation on 

Hp ermil which implies that (n;U,t) <-^> (1 — n; —U,t). 
The present approach also gives a reasonable picture for 
U = 0. In this case, the ground state energy of the bath 
is exactly half of the energy of the half-filled Fermi-sea, 
because the spinless fermions do not account for the elec- 
tron spin. In principle, the remaining half could be re- 
trieved from Hp au u, if each nearest-neighbor bond were 
in the singlet state. It clearly points at the singlet charac- 
ter of the free-electron ground state. Also, the J is linear 
in t (unlike the superexchange J), and gives a measure 
of the exchange energy of the Fermi-sea. It appears to 
us that this representation may be able to describe the 
large and small U limits in a unified way. We will present 
more studies on Hubbard model in a future work. 



IV. \U\ = oo: NAGAOKA FERROMAGNETISM 

We now address the case of infinite U Hubbard model. 
It turns out that taking U — * ±oo limit is straightforward 
in the new representation. Since U is the chemical po- 
tential for spinless fermions, the infinite \U\ completely 
suppresses the spontaneous fluctuations in the number 
of spinless fermions. That is, in the limit \U\ — » oo, 
the grand canonical bath of spinless fermions becomes 
canonical. For positive-infinite U, the Hubbard model 
on square lattice is known to have ferromagnetic ground 
state close to half- filling (Nagaoka ferromagnctism 25 ), 
which eventually becomes unstable around n e ~ 1 ± 8 C 
(early estimate of the critical (hole) concentration is 
5 C ~ 0.49 34 ; later works suggest S c ~ 0.25 35 ). There have 
been many studies on the existence and stability of the 
Nagaoka ferromagnetism on different types of lattices. 
Presently, we discuss this phenomenon in our representa- 
tion, for positive as well as negative infinite-?/ Hubbard 
model on bipartite lattices. The \U\ = oo Hamiltonian 
in this case can be written as: 



-ffoo — —t 



[l + <xr • TR+s) 



*R^R+5 



&R+<5 a 



R 



(17) 



where the hopping of the spinless fermions is exclusively 
determined by the exchange operator on each nearest- 
neighbor bond, in a way, similar to (but not same as) 
the Anderson-Hasegawa double-exchange 36 in a ferro- 
magnetic metal. Although i/oo is the same for both 
signs of infinite-?/ , the two limits correspond to differ- 
ent states of the spinless fermion bath. For U = +oo, 
the bath is fully- filled (n = 1), while it is completely 
empty (n = 0) for U = — oo. Both of these cases are 
dead, as there are no fermions which can gain energy 
by hopping. Also, there is 2 L -fold degeneracy due to L 
independent Pauli objects. In terms of electrons, the for- 
mer (U = +oo) corresponds to the localized electrons at 
half-filling, while in the latter case, an arbitrary site in 
a configuration is cither empty or doubly occupied. In 
order to discuss interesting physical possibilities, we 'ex- 
ternally' populate (or depopulate) the fermion bath, by 
employing the chemical potential, (i, for electrons. That 
is, we consider /Coo = Hoc— fi^N^+Ni), where N^+Ni = 

Er^II + (1 - «r«r)41 + E ReB [l + (1 - blb R )r^}. 
Below, we shall discuss the existence and stability of the 
ferromagnetic ground state in /Coo. We will see that the 
present approach is considerably simpler as compared to 
the usual Gutzwiller projection based methods. 

First we discuss Hoc on just two sites. The exchange 
operator on a bond has eigenvalue +1 for the triplet and 
— 1 for the singlet state of the two-level objects. Al- 
though the complete eigenstates of Hoq in the two cases 
are degenerate, they physically presents two distinct pos- 
sibilities in the ground state. For a fully polarized 
triplet |— , — ), the ground state of the two-site problem 
is: (|l,0) + |0,l))<g> |-,-) = | 1,0) + 10,1) (ignoring the 
normalization factors). This is a trivial 'ferromagnetic- 
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metallic' state with J, magnetization. Similarly for | + , +), 
it is | t; Tl) + I t!>t)j which is also ferromagnetic (but 
with respect to the fully- filled electron configuration). 
For the singlet state, (|+, — ) — |— , +)), the ground state 
is: (|1,0) - |0,1)) ® (|+,-) - |-,+)) = | t,0) + |0,T 
) + I T4j I) + I Jo Tl}) which is like a 'non-magnetic' metal- 
lic state. Also, in the singlet case, n-f = nj. = 1/2 and 
S = 0. This simple analysis suggests two qualitatively 
different possibilities in the ground state of H^. 

Now we study JC^ on the full lattice. We explicitly 
consider the case of U = — oo (that is, gradually fill the 
empty bath) . The results for U = +oo can be exactly in- 
ferred from the particle-hole transformation on the spin- 
less fermions, which states that the ground state for +ve 
U for a given density n, is same as that of —we U for 1— n, 
for a fixed t and m [where m = (J2a °^ +S_b tZ )I^\- To 
establish the existence of ferromagnetic-metallic ground 
state for U = — oo, consider one spinless fermion in the 
empty bath. The (only) energy, that it gains by hop- 
ping, can be maximized by having maximum and uniform 
hopping through each bond. While the fully-polarized 
'ferromagnetic' state of the Pauli subsystem fulfills this 
requirement, the same can not be achieved in any other 
state, for instance, a configuration of bond-singlets or any 
deviation with respect to the fully polarized state. The 
fully polarized Pauli subsystem therefore emerges as the 
only choice which helps the fermion in gaining maximum 
kinetic energy. The same will be true for a finite density 
of the spinless fermions, until the 'ferromagnetic' state of 
the Pauli subsystem becomes unstable. In terms of elec- 
trons, the ground state of a finite n fermion bath together 
with '— ' polarized Pauli subsystem, corresponds to a J, 
polarized ferromagnetic metal with the electron density, 
n e = n. Also, for the '+' polarized Pauli subsystem, the 
electronic ground state is ferromagnetic, but from the 
fully-filled side of the electron density (n e = 2 — n). This 
analysis further implies that the ground state ofU = +oo 
Hubbard model is also ferromagnetic-metallic, as discov- 
ered by Nagaoka, around half-filling (n e = 1 ± n, for '+' 
or ' — ' polarization respectively). Below we discuss the 
instability of the ferromagnetic state. It is generally a 
hard problem. But in our framework, it actually becomes 
a much simpler problem. 

As it appears from the structure of Hoc, the spinless 
fermions are happy to live with ferromagnetic Pauli sub- 
system for any density of electrons. But the Pauli sub- 
system itself is not comfortable with the change in n e , be- 
cause \i acts as an 'external' field for the Pauli operators. 
This induces 'flips' in the fully-polarized Pauli subsystem 
for sufficiently populated fermion bath, thereby causing 
the instability of the Nagaoka state. To work this out, 
we write the Pauli operators in terms of d and (3, the two 
hard-core bosons, such that: a z = — 1 + 2dra, a + = cr 
and t z = — 1 + 2/3' 0, r + = fy , The creation of a hard- 
core boson amounts to flipping |— ) on the respective site 
to |+), like a 'magnon'. Rewriting /C^ in terms of a and 
separates the ferromagnetic part from the terms which 
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FIG. 3: Metallic ferromagnetism in the infinite- U Hubbard 
model. The shaded region between n e = 0.5 and 1.5 is the 
Nagaoka ferromagnetic state for U = +oo. Note the com- 
plementary existence (in the sense of electron filling) of the 
Nagaoka ferromagnetism for U — +oo and — oo. 



cause flips, that is: /Coo = H FM + , where 

REA 5 

TLEA RGB 

and 

H fUps = t J2Y, (aiPll+S + ^R+5«r) X ( 19 ) 
RSA 5 

[(Q4a R + P^+Jn.+s) - (a R /3 R+ «5 + /3 r+ 5"r) 

~ 2 M ( X ~ «R«R)«R d R 
Re A 

-2H J2 (1 - &rM/3r/3r 
reb 

In the ferromagnetic state, H f lips is zero as all the two- 
level objects are in the same state, and only H FM ac- 
counts for the ground state for different n e , where n e = n. 

Instability of the ferromagnetic state can be deter- 
mined by calculating the dispersion of exactly one 'flip', 
in the ground state of H FM for different values of /i. 
That is, we replace the fermion operator terms in H fi ips 
by their expectation values in the ground state of H FM , 
and calculate the one- flip dispersion. That value of n e , 
for which the minimum of this dispersion becomes nega- 
tive, marks the instability for the Nagaoka state. By fol- 
lowing the above prescription, we get the following equa- 
tions for the spinless femion density, n, and the one-flip 
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dispersion, e^. 

n = 7 E [ e (~^l7k| - M) + e(2i| 7 k| - M)] (20) 



k 



(21) 



Y, bkl [e(-zi| 7k | - m) - e(zt|7k| - m)] 



e± = ^ X (l±|7k|)-2/x(l-n) 



(22) 



Here, % is the ferromagnetic ground state energy per 
nearest- neighbor bond (in units of — t), and tx is always 
positive (for any sign of t). Since ztx(l ± |7k|) > and 
(1 — n) is also positive, therefore, ej^ is certainly non- 
negative for fi < 0, which implies a stable ferromagnetic 
ground state. For fi = + however, becomes negative 
at k = 0. Thus, the ferromagnetic ground state becomes 
unstable to the flips in the Pauli subsystem, and there 
occurs a first order transition in the ground state. The 
stable ferromagnetic-metallic ground state (/x < 0) cor- 
responds to n e = n < 1/2, and the instability exactly 
occurs at n e = 1/2 (zx = 0). The same analysis for '+' po- 
larized case gives the stable Nagaoka state for n e = 2 — n, 
where < n < 1/2. Using the particle-hole transforma- 
tion for spinless fermions, we further deduce that the 
Nagaoka state for U = +oo exists for 1/2 < n e < 1 (hole 
doping) and for 1 < n e < 3/2 (electron doping). These 
results are presented in Fig. 3. 

This simple analysis of a difficult problem is very en- 
couraging. It gives the essential physics, that is the fer- 
romagnetic ground state and its instability, in a very 
straightforward manner. We could show the existence of 
Nagaoka state for both U = +oo as well as — oo Hubbard 
model. Of course, we would like to investigate (more ac- 
curately) whether the critical n is exactly 0.5 or slightly 
different. But n = 0.5 appears to be a special fermion 
density for bipartite lattices (maximum ferromagnetic 
exchange; also, corresponds to a non-magnetic metal- 
lic ground state). Further investigations of the Nagaoka 
ferromagnetism for finite U case, and on non-bipartite 
lattices will be discussed elsewhere. Finally, an impor- 
tant point to note in the present discussion is that the 
ferromagnetism in Hubbard model is strictly related to 
the invariance of the spinless fermions under a global 
gauge transformation, while the antifcrromagnctic/non- 
magnetic states break this gauge invariance. This is a 
novel point of view for the strongly correlated electrons. 



V. SUMMARY 



as the chemical potential for spinless fermions. 
On a bipartite lattice, the Anderson's superex- 
change for large and positive (as well as negative) 
U is shown to correctly emerges from the spinless 
fermion bath, within a simple decoupling scheme. 

• A Hamiltonian is derived for studying Nagaoka fer- 
romagnetism in \U\ — oo Hubbard model. The ex- 
istence of the ferromagnetic-metallic state and its 
stability is discussed. It is pointed out that the fer- 
romagnetism in Hubbard model respects the global 
gauge invariance for the spinless fermions, while the 
antiferromagnetism necessarily violates it. 

• The meaningful analysis of the two important cases 
of the Hubbard model is very encouraging for the 
new representation, and demonstrates its useful- 
ness for studying strongly correlated electrons. 

• A representation for bosons is discussed in ap- 
pendix B. The Jordan- Wigner transformation is 
derived in appendix C by introducing a unitary 
transformation with a Majorana fermion as gauge. 
Finally, the notion of a soft-core fermion is intro- 
duced in appendix D. 
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APPENDIX A: GENERALIZATION TO 
MULTICOMPONENT FERMIONS 

Let fj be i th component of an iV-component fermion 
operator. Its Hilbcrt space is 2 N dimensional. In terms of 
a spinless fermion and N — 1 two-level objects, both the 
equality of the Hilbert space dimensions and the anti- 
commutativity can be attained consistently for the N- 
component fermion, by defining a representation in the 
following way. 



fl 
fl 

f 

JN 



i=i 



■J at for i = 2, • • 



■ 7 1 N ~ 1 



(Al) 

,7V - 1 (A2) 
(A3) 



i=l 



• A new canonical and invertible representation for 
the electrons, in terms of a spinless fermion and the 
Pauli operators, is presented. It is further general- 
ized for the TV-component fermions in appendix A. 

• The Hubbard interaction in this representation acts 



The operators fj (for i = 1,N — 1) anticommute with 
each other because anticommutes with erf, and fj^ 
anticommutes with the rest of fj operators, because tp 
anticommutes with <j) (in addition to the anticommuta- 
tion of (T 4 + and erf). The inverse of this representation 
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[Eqs. (A1-A3)] is given as: 
For 



i = l,N-2 

N-l 

f} J] {2n 3 - 1) (/ v + /, 
j=i+i 



N 



N-l — /jV- 


i 


(Pn + In) 


« f = II 


1 


-nf=T 1 (2« i -i) 




2 


In 


1 


+ n£ 1 (2ni-l) 


2 



(A4) 
(A5) 



(A6) 



where hi = fjfi. The similarity between Eqs. (A4- 
A5) and the Jordan- Wigner transformation is notable, 
with the exception that the Jordan- Wigner does not 

have an extra (fj^ + f^) operator [and no equivalent 

of Eq. (A6)]. Moreover. Eqs. (Al) and (A2) are sim- 
ilar to the inverse Jordan- Wigner transformation [with- 
out a <fi and Eq. (A3)]. One can generate many equivalent 
representations for the TV-component fermion by making 
suitable unitary rotations 37 . 



APPENDIX B: 



REPRESENTATION FOR 
BOSONS 



Let and a be the creation and annihilation operators 
of a canonical boson, acting on the Fock states, {\n), n = 
0,1..., oo}. The operator can be written as: 



Here, a are the Pauli operators defined in the Hilbcrt 
space of |e) and |o), such that er+ = |o)(e|. Operators W 
and b are the new canonical bosons, acting in the Fock 
space spanned by {|m)}, which commute with the Pauli 
operators. It is easy to check that Eq. (B5) satisfies the 
bosonic commutation relations. The harmonic oscillator 
in terms of the original bosons takes the following form 
in the new representation. 



a)a+ i = 2 



6*6- 



1 



(B6) 



The new form looks like a two-level atom in a quantized 
single-mode radiation field (without atom- field coupling) , 
which conversely suggests that an atom-field problem can 
be exactly mapped to an equivalent field-only problem (a 
kind of 'unification' of matter and radiation-field into a 
new radiation-field). Most direct and simple case of this 
unification, present in Eq. (B6) itself, is when the fre- 
quency of radiation is twice the energy-difference between 
the two atomic- levels [see right-hand side of Eq. (B6)]. 
The presence of an atom-field interaction will further in- 
troduce non-trivial terms in the unified radiation field. 
The application of this representation to an atom-field 
interaction problem will be discussed elsewhere 38 . 

Let us write down the inverse of this representation. 
The Pauli operators, in terms of a and a^, are given as: 



cos7ra^a 
1 



i-x 

2 



1 



(B7) 
(B8) 



^2Vn~TT \n + l)(n\ (Bl) 

n=0 
oo 

\/2m+ 1 |2m+ l)(2m| (B2) 

m=0 

oo 

+ ^ y/2(m+l) |2(m + l))(2m+l| 



where n = 2m in the first summation in Eq. (B2), and 
n = 2m+l for the second. This even-odd structure of the 
Fock states allows us to write a composite representation 
of the boson operators by defining the following mapping. 



1 2m) := \m) ® \e) 
|2m+ 1) := |m) <g> \o) 



(B3) 
(B4) 



The kets |e) and |o) stand respectively for the even and 
odd values of the original Fock state quantum number n, 
and the states {|m)} define a 'new' Fock space. Substi- 
tuting this mapping in Eq. (B2) results in the following 
representation of a* . 



2tfb+ 1 cr H 



V2 sv 



(B5) 



Note that (1 ± x)/2 are the projection operators of the 
odd and even states respectively, in the original Fock 
space. Therefore, the meaning of Eq. (B8) is clear that 
it connects an even \n) to an odd |n) and annihilates the 
odd states, consistent with the definition of a + . Also, 
the Pauli operators in the above representation satisfy 
the necessary algebra. Interestingly, Eqs. (B7) and (B8) 
also give us an unconstrained new bosonic representation 
for spin- 1/2 operators. If we restrict the Fock space to 
|0) and |1) (hard-core constraint), then we retrieve the 
usual hard-core boson representation. Now to the boson 
which, in terms of a and a) , is represented as: 



^ = 4= 

V2 
1 

V2 



i-x 



1 



i-x 



— - — a'a' 

2 VdKTT 



(B9) 



It is consistent with the fact that the number of 6-bosons 
changes by 1 when the number of a-bosons changes by 2. 
Besides, it satisfies the bosonic commutations, and com- 
mutes with a z and , as in Eqs. (B7) and (B8). Hence, 
a canonical and invcrtible representation for bosons. 
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APPENDIX C: UNITARY TRANSFORMATIONS 
WITH 'MAJORANA' GAUGE 

We now introduce a new kind of gauge transformation, 
and show how it generates the Jordan- Wigner represen- 
tation (JW) for spin-1/2 operators. Usually, we know the 
JW as an intuitive mapping that was not derived. Here, 
we precisely do that: derive it. Consider er + in terms 
of the electron operators [as in Eq. (5)]. It is given as: 
o" + = where , = ft + f\ is a Majorana fermion 



corresponding to | electron. Since 



1, therefore, 



±e l 2 (j). , and a + can be rewritten as: 



(CI) 



where the [ electron Majorana fermion appears in the 
gauge. This observation motivates us to define a suitable 
'Majorana' gauge-transformation such that the phase 
factor in Eq. (CI) is completely absorbed into /|, thereby 

transforming o~ + into ft. We easily define such a trans- 
formation as: 



-£-(1 + ^)7^ 



U(<f)^) = exp 

= (1 - n T ) - 0_ L n t 



(C2) 



where n^. = /|/-p and the expression in the second line 
has been derived from the first by using the properties: 

(^IT 1 ) = ~~T^ an d "t = Clearly, U(4>±) is both 
unitary as well as Hermitian. It is a simple extension of 
the usual gauge transformation, achieved by allowing a 
Majorana fermion to be the phase. One can show now 
that W^^cr+Ui^) = /j, while a z remains the same, 
that is 2rt| — 1. This is nice but unusual, because under 
a 'normal' unitary transformation, a canonical fermion 
transforms into another canonical fermion. But here, 
under the Majorana gauge transformation, a canonical 
fermion becomes a Pauli operator (hard-core boson) and 
vice-versa. This change in the canonical nature happens 
due the fact that the gauge itself is a real fermion, not just 
a number. Next, we discuss the consequence of applying 
the Majorana gauge transformation [as in Eq. (C2)] on a 
collection of spin-1/2 objects. 

Let there be L number of spin-1/2 objects described 
in terms of the Pauli operators, cr ; + , where the integer 



label I goes from 1 to L 
?t 



In order to transform every a. 



into a corresponding fL, we must apply W(<j^) for each I. 

Since U{<f>i,)U{4> m ,) ^ U{(f) m ^)U{(f) l ^), it is important to 
keep track of the order in which to make the local gauge 
transformations. We can do this by defining an ordered 
unitary operator, U, as follows. 



u=u(<j> Ll ) ■■■ u{4> 2i )u(<j> n ) 



(C3) 



In principle, we could define Li in any order, but we just 
choose to work with the above arrangement. Now, we 



TABLE III: Algebra of the soft-core fermions, and its com- 
parison with that of the hard-core bosons. Here, [ , ] denotes 
a commutator, and { , } denotes an anticommutator. 



uUl I vUI c .T cl IlllUIlo 


lldl LI L/UI U JJUSUIIn 


\y- y+1 _ i 

[Zj , Zj J — 1 


In' n + \ - 1 

X" i u J — 1 


yz _ ly+ y- \ 

Zj — \Zj , Zj J 


a * _ r_+ _-] 


\y z y^] — \-9Y ± 


\„* ±1 _ +2 ± 
\ U . (J — ZL Zi U 


Hilbert Space: oo dimensional 


2 dimensional 


Basis: {|n} : n — 0, 1, . . . } 


{|+>,|-» 


S z |n) = (2n + l)|n) 


<r*|±> = ±|±> 


E+|n) = y/n + l\n + l) 


^ + l-) = l+) 


E~|0) = 


a~\-) = 



apply hi one by one on each ay , and get the following 
answer, which is the JW transformation. 



WV+W = f\ 



T 
i-i 



V I 



L 



■hi 



(C4) 
(C5) 



We have thus derived the JW representation for spin- 
1/2 operators by applying a unitary transformation (with 
Majorana gauge) on the Majorana representation of the 
same. The origin of the string of (1 — 2h^) in the JW 
representation also becomes clear from the nature of U. 
For a ; + , the operators U((j) L ,) to U((j> l+1 ^) do nothing, 
while UQ>\) gives fL. The rest which follow afterwards, 
in the given order, contribute a factor of (1 — 2?/L) each, 

because every U((j> m ^) (V m < I) that commutes past f^ 
becomes U{-$ mi ), and U(<p ml )U{-<p ml ) = 1 - 2n mV 



APPENDIX D: SOFT-CORE FERMION 

In this appendix, we introduce a new object, called 
soft-core fermion, which is an anticommuting counter- 
analog of the hard-core boson (the Pauli operators). We 
define the raising operator S + of a soft-core fermion as: 

£+ = tf<j> (Dl) 

where w is the creation operator of a canonical boson, 
and <f> as usual denotes a Majorana fermion such that 
<j> = & + c, where c and are the operators of a canoni- 
cal fermion. This definition of a soft-core fermion is mo- 
tivated by the Majorana representation [Eq. (5)] of a + , 
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the hard-core boson. As the name suggests, a soft-core 
fermion anticommutes with other fermions and soft-core 
fermions, while behaving as a boson for itself. This lat- 
ter property implies no local constraint on the number 
of such particles (opposite of the hard-core constraint), 
hence the name soft-core. Also, it commutes with other 
bosons and hard-core bosons. The algebra of the soft- 
core fermions is presented in Table III, and compared 
with the hard-core bosons. 

Interestingly, we can derive a JW type representation 
for the soft-core fermions, in exactly the same way as we 
did for the Pauli operators in the previous appendix. We 
define a unitary transformation: 



u{4>) 



exp 



-4(1 



(D2) 



ri + (-)*- 






2 


-4> 


2 



such that U^^+U^) = &t. Here, h = b% and (-)" = 
e lnn = cos (irh) . If n were the number operator of a 
fermion, then Eq. (D2) correctly reduces to Eq. (C2). 
Next, we define an ordered unitary operator hi for L such 
soft-core fermions, U = U{cf)i)U(<j)2) ■ ■ -U^l)- Under U, 
the soft-core fermions transform as: 



(-) 
V I = 2 



(D3) 
(D4) 



which is same as the JW transformation, but in terms 
of the bosons. We have thus presented the mathemati- 
cal structure of soft-core fermions, whose actual physical 
usefulness is not obvious at present, but might become 
so in future. 
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